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ABSTRACT 

The prompt emission spectra of gamma-ray bursts (GRBs) usually have a dominant component 
that is well described by a phenomenological "Band" function. The physical origin of this spectral 
component is debated. Although the traditional interpretation is synchrotron radiation of non-thermal 
electrons accelerated in internal shocks or magnetic dissipation regions, a growing trend in the com- 
munity is to interpret this component as modified thermal emission from a dissipative photosphere 
of a GRB fireball. We analyze the time dependent spectrum of GRB 110721A detected by Fermi 
GBM and LAT, and pay special attention to the rapid evolution of the peak energy E p . We define a 
"death line" of thermally-dominated dissipative photospheric emission in the E p — L plane, and show 
that E p of GRB 110721A at the earliest epoch has a very high E p ~ 15 MeV that is beyond the 
"death line" . Together with the finding that an additional "shoulder" component exists in this burst 
that is consistent with a photospheric origin, we suggest that at least for some bursts, the "Band" 
component is not from a dissipative photosphere, but must invoke a non-thermal origin (e.g. syn- 
chrotron or inverse Compton) in the optically thin region of a GRB outflow. We also suggest that the 
rapid "hard-to-soft" spectral evolution is consistent with the quick discharge of magnetic energy in a 
magnetically-dominated outflow in the optically thin region. 

Subject headings: gamma rays:bursts — gamma rays: observations — gamma rays: theory — plasmas — 
radiation mechanism: non-thermal — radiation mechanism: thermal 



1. INTRODUCTION 

The prompt emission spectrum of a gamma-ray burst 
(GRB) is usually well described by a phenomenologi- 
cal function known as the "Band" function (Band et al. 
1993). This model, which is essentially a broken power 
law function with a smooth (exponential) transition, was 
traditionally invoked to model spectra of GRBs detected 
by BATSE on board the Compton Gamma- Ray Observa- 
tory. The function is found successful to describe most 
GRB spectra detected by later missions as long as the 
spectral band is wide enough (e.g. Abdo et al. 2009; 
Zhang et al. 2011). 

The physical origin of this phenomenological Band 
function is not identified. The traditional model is syn- 
chrotron emission of non-thermal electrons in an opti- 
cally thin region, e.g. internal shocks or internal mag- 
netic dissipation regions (Meszaros et al. 1994; Tavani 
1996; Daigne & Mochkovitch 1998; Lloyd & Petrosian 
2000; Bosnjak et al. 2009; Zhang & Yan 2011). Alter- 
natively, a matter-dominated outflow (fireball) can have 
a bright photosphere (Paczynski 1986; Goodman 1986; 
Meszaros & Rees 2000; Meszaros et al. 2002), which 
may be enhanced by kinetic or magnetic dissipation pro- 
cesses near the photosphere (Thompson 1994; Rees & 
Meszaros 2005; Pe'er et al. 2006; Giannios 2008; Be- 
loborodov 2010; Lazzati & Begelman 2010; Ioka 2010). It 
has been argued that due to geometrical and / or physical 
broadening, this quasi-thermal component may be mod- 
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ified to mimic a Band function (e.g. Beloborodov 2010; 
Lazzati & Begelman 2010; Pe'er & Ryde 2011; Lundman 
et al. 2012). The scalings of a non-dissipative photo- 
sphere model are argued to be able to interpret various 
empirical correlations (Fan et al. 2012). 

Within the framework of the standard fireball-shock 
model, a GRB prompt emission spectrum is expected 
to be the superposition of a quasi-thermal photosphere 
emission component and a non-thermal component in 
the optically-thin internal shock region (Meszaros & Rees 
2000; Zhang & Meszaros 2002; Toma et al. 2011; Pe'er et 
al. 2012). Such a superposition effect has been claimed 
in the BATSE data archive (e.g. Ryde 2005; Ryde & 
Pe'er 2009), and was confirmed more robustly recently 
with the Fermi data (e.g. Ryde et al. 2010; Zhang et al. 
2011; Guiriec et al. 2011; Axelsson et al. 2012). On the 
other hand, most GRB spectra (e.g. GRB 080916C) are 
still well described by one single Band component (Abdo 
et al. 2009; Zhang et al. 2011). This sharpens the de- 
bate regarding the origin of the Band function. For GRB 
080916C, the non-detection of a thermal component led 
to the suggestion of a Poynting-flux-dominated outflow 
(Zhang & Pe'er 2009; see also Daigne & Mochkovitch 
2002; Zhang & Meszaros 2002). Alternatively, some au- 
thors attempted to interpret the entire Band function 
as emission from a dissipative photosphere (e.g. Be- 
loborodov 2010; Vurm et al. 2011; Giannios 2008; Ioka 
2010). These two interpretations invoke distinct assump- 
tions regarding the composition of the GRB jets. Find- 
ing observational clues to differentiate between them is 
therefore essential to unveil the physics of GRB central 
engine, jet composition and energy dissipation mecha- 
nisms, which are poorly constrained (e.g. Zhang 2011). 

Here we show that the time-resolved spectral informa- 
tion of GRB 110721A holds the key to address this open 
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question. 

2. "DEATH LINE" OF GRB BARYONIC PHOTOSPHERE 
EMISSION IN THE E P — L PLANE 

For a hot fireball with total wind luminosity L w 
launched from an initial fireball radius Ro, the initial 
temperature is 

T ~ (LjAirRlca) 1 / 4 ~ 1.4 x 10 10 K L^R^ 2 , (1) 

where c is speed of light, and a = 7.56 x 
10~ 15 erg cm -3 K -4 is the Stefan-Boltzmann energy 
density constant. The observed photosphere tempera- 
ture T p h can be as high as To if the fireball is clean enough 
so that the photosphere radius R p h does not exceed the 
fireball coasting radius R c , but is lower than T for fire- 
balls with a heavier baryon loading when R p h > R c - 
More specifically, one has (Meszaros & Rees 2000) 

^ = { (^)~ 2/3 =fe) 8/8 ' r,<v*,Rph>Rc, 
T ° I 1, V > r)*,Rph < 

(2) 

where rj = L w / Mc 2 is the dimensionless entropy of the 
fireball, and 

is the critical value of rj. The photosphere luminosity 
is L ph ~ tt(1 /T) 2 R 2 aT 4 oc (R/T) 2 T 4 oc R 2 T 2 T' 4 . For 
77 > ?/*, since T oc R and T' oc R^ 1 , one has oc R° ~ 
const. For 7/ < 77*, since T oc T oc i?~ 2 / 3 , T^ oc 
R 2 R~ 8 / 3 oc i?~ 2 / 3 . So the photosphere luminosity is 

r f /B \- 2 / 3 / \ 8 / 3 

^ = J (iff) =(i) . v<v*,R P h>Rc, 
Lw I 1, 17 > < -Rc, 

(4) 

If a GRB spectrum is dominated by the photosphere 
emission, for a certain observed isotropic 7-ray luminos- 
ity L — L p h, the spectral peak energy E p should not 
exceed (kT , where ( is a factor to denote the vF v peak 
of the photosphere spectrum. Therefore a baryonic pho- 
tosphere emission has a "death line" defined by 

E p < (kT ~ 1.2 MeV (L 1 ^ 4 R^ 2 . (5) 

The factor £ is subject to the shape of the spectrum. For 
a strict blackbody, C ~ 3.92. For a relativistic outflow, 
the shape of blackbody is modified to the form (see also 
Li & Sari 2008 and references therein) 

v 2 f°° dx 

c Jhv/kT e - 1 

The vF v spectrum has a maximum value at 

C - 2.82. (7) 

The death line ([5]) is derived for a non-dissipative pho- 
tosphere. For a thermally-dominated jet, dissipation 
via internal shocks or neutron collisional heating near 
the photosphere would compensate adiabatic cooling, so 
that the photosphere temperature would be maintained 



close to but never exceed the maximum temperature (Be- 
loborodov 2012), so eq.© applies to a broad categories 
of dissipative photosphere models as well. A possible 
exception would be the dissipative photosphere model 
that invokes continued magnetic dissipation at extended 
radii above the photosphere (Drenkhahn & Spruit 2002), 
which allows E p to be above the death line (J5J) if the bulk 
Lorentz factor is large enough and dissipation at high op- 
tical depth is suppressed (Giannios 2012). However, in 
order to thermalize the jet, the required Lorentz factor 
is very low (Vurm et al. 2012), making it essentially im- 
possible to exceed the death line. Also this model may 
not interpret the GRB 110721A phenomenology, includ- 
ing the existence of the shoulder thermal component and 
the rapid spectral softening during the rising phase of 
bolomctric luminosity (see Sect. 3). 

3. GRB 110721A 

GRB 1 10721 A was jointly detected by the Gamma- Ray 
Burst Monitor (GBM; Meegan et al. 2009) and the Large 
Area Telescope (LAT; Atwood et al. 2009) onboard the 
Fermi Gamma-Ray Telescope (Axelsson et al. 2012 and 
references therein). A candidate optical counterpart was 
reported (Greiner et al. 2012). Assuming that the as- 
sociation is real, Berger (2012) suggested two possible 
redshifts z = 0.382 or z — 3.512, with the former one 
preferred. 

The time-dependent spectral evolution of GRB 
110721A was reported by Axelsson et al. (2012). Two 
noticeable features of the burst are: 1. A thermal com- 
ponent is identified to be superposed on the Band com- 
ponent in both the time integrated spectrum and time- 
resolved spectra. The temperature of this "shoulder" 
thermal component evolves with time as a broken power 
law, which is consistent with the expectations of the pho- 
tosphere model (Ryde & Pe'er 2009). 2. The Band com- 
ponent displays an extremely rapid spectral evolution. 
The E p at the earliest epoch reaches a record-breaking 
value ~ 15 MeV. 

We have independently processed the Fermi data of 
GRB 110721 A. We performed a joint spectral analysis 
using the data from the Nal 6, 7 and BGO 1 detec- 
tors on GBM, as well as the LAT dat£0. For GBM, we 
used the TTE event data containing individual photons 
with time and energy tags. Background rates are esti- 
mated by fitting the light curve before and after the burst 
using a one-order background polynomial model. We 
pretreated the LAT data using the LAT ScienceTools- 
v9r27pl package and the P7TRANSIENT_V6 response 
function (detailed informatin for the LAT GRB Analy- 
sis are available in the NASA Fermi web sit^J). In the 
diffuse response calculation, a three-component model 
is used: GRB 100721 A with a PowerLaw2 spectrum, 
the Galactic diffuse model of gal_2yearp7v6_v0.fits, and 
the extragalactic diffuse power law model. We then ex- 
tracted the background-subtracted light curves from the 
GBM and LAT and demonstrated them in Figure [T] It 
is obvious that higher energy photons arrive earlier than 
the lower energy photons. 

4 http:/ /fcrmi. gsfc.nasa.gov/cgi-bin/ssc 

/LAT/LATDataQuery.cgi 

http : / / fermi .gsfc. nasa.gov /ssc /data/ analysis /scitools 
/lat_grb_analysis . html 
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Fig. 1. — The light curves of GRB 110721A in different energy 
bands: solid: Nal; short dash: BGO; long dash: LAT above 100 
MeV. 

In order to better understand the cause of such a be- 
havior, we carry out a detailed time-dependent spectral 
analysis. We adjust the size of the time bins so that 
each time bin contains enough photons to perform a sta- 
tistically significant spectral analysis. We applied the 
software package RMFIT (version 3.3pr7) to carry out 
the analyses. We confirm the conclusion of Axelsson et 
al. (2012) that adding a shoulder thermal component 
can improve the fits to the data significantly. Since the 
evolution of the thermal component has been presented 
in Axelsson et al. (2012), in this paper we focus on the 
Band component only in order to study its physical ori- 
gin. We find that the Band model usually gives a rea- 
sonable fit to the data, with the reduced \ 2 i n the range 
~ (0.9 — 1.1). Figure [2] gives an example of the Band 
model fitting the data in the time interval of [-0.512, 
0.064] s. This earliest epoch indeed shows an extremely 
high E p value ~ 19.6 ± 4.5 MeV, which is consistent 
with E p = 15 ± 1.7 MeV reported by Axelsson et al. 
(2012). The larger error in our fit may be because the 
Fermi team has included the extra LLE (LAT Low En- 
ergy) data, which is currently unavailble to the public. 
We also present the best fit vF v model curves in different 
time intervals in Fig|3] 

Figure [4] displays the rest- frame E p — L plot for time- 
resolved spectra of Fermi GRBs with well measured red- 
shifts. The data are a sub-sample of Fig. 9 of Lu et al. 
(2012), for which only the parameters during the rising 
phase of GRB pulses are adopted. This is because the 
decaying phase may be controlled by the high-latitude 
curvature effect, which does not directly reveal radiation 
physics. Two "death lines" (eq.[5]) are drawn, which 
correspond to Rq ~ 10 7 cm (solid, typical value) and 
i?o ~ 3 x 10 6 cm (dashed, an extreme value to allow 
highest death line possible). GRB 110721A at the ear- 
liest epoch [-0.512. 0.064]s is plotted for two candidate 
rcdshifts. It is clearly seen from the figure that for both 
redshifts, the points are way above the death lines. This 
rules out a wide range of dissipative photosphere models 
at least for this earliest epoch. 

One may make one step further. Axelsson et al. (2012) 
showed that E p evolution of the Band component can be 
fit as a power-law decay with time, suggesting that this 
component might have a same physical origin during the 
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Fig. 2. — Data and Band function model curve of GRB 110721A 
in the time interval of [-0.512, 0.064]s. The solid line is the 
best Band function fit with parameters a = —0.944 ± 0.046, 
= -4.543 ± 0.887, and peak energy E p = (1.956 ± 0.422) X 10 4 
keV, which are consistent with Axelsson et al. (2012). 
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Fig. 3. — The best-fit vF v model spectra for the time resolved 
data in different time bins (as marked in the legend). 




Fig. 4. — The rest-frame peak energy E p (l + z) is plotted against 
the observed isotropic 7-ray luminosity L for the time resolved 
spectra of Fermi GRBs with redshift measurements. The black 
data points are from Fig. 9 of Lu et al. (2012), but only the 
spectra during the rising phase of GRB pulses are taken (see text 
for explanation). The blue points are the time-resolved spectra of 
GRB 110721A for the two candidate redshifts. The two stars are 
for the first epoch, which are both well beyond the death lines. 
Two death lines are plotted, which correspond to Rq = 10 7 cm 
(solid) and 3 X 10 e cm (dashed), respectively. 
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burst. Figure|3]displays evolution of the spectral shape of 
the Band component, which shows a similar a value with 
a gradually shallowing (3 value. This suggests that the 
Band component should share the same physical origin 
in different epochs: it is not from the photosphere, but 
is likely from an optically thin region where non-thermal 
particles are accelerated. 



4. SUMMARY AND DISCUSSION 

We have shown that the Band function component of 
GRB 110721A is beyond the "death line" of thermally- 
dominated dissipative photosphere models in the E p — L 
plane. Together with the fact that an additional shoulder 
thermal component is consistent with the photosphere 
model, we reach the conclusion that the so-called "Band" 
component is not of a photospheric origin at least for 
this burst, and is formed via non-thermal dissipation pro- 
cesses in the optically thin regions. Veres et al. (2012) 
interpreted this emission as synchrotron emission from a 
magnetically-dominated jet. 

Such a finding has profound implications in under- 
standing the origin of other Band function spectra of 
GRBs. Zhang et al. (2011) identified three elemen- 
tal spectral components through a detailed time-resolved 
spectral analysis of 17 GRBs co-detected with Fermi 
GBM and LAT. They found that there are two types 
of GRBs. GRB 080916C's spectra remain the "Band" 
shape even though the time interval for the spectral anal- 
ysis progressively reduces, reaching ~ 1 s in the rest 
frame. GRB 090902B, on the other hand, showed a clear 
"narrowing" feature as time bin reduces, and spectrum is 
much narrower at ~ rest-frame 1 s. The time integrated 
spectrum of this burst is also narrower than other Band 
GRBs. This burst's "Band" component can be indeed 
de-composed as superposed photosphere emission (Ryde 
et al. 2010; Zhang et al. 2011; Pe'er et al. 2012; Mizuta 
et al. 2011). However, such bursts are not common. 
Most bursts are similar to GRB 080916C. 

The spectral parameters of GRB 110721A are similar 
to those of GRB 080916C and many other GRBs. The 
conclusion that the Band component of GRB 110721A 
originates from the optically thin region also supports the 
suggestion that most Band components are non-thermal 
(synchrotron or SSC) emission in the optically thin re- 
gion. Available data seem to suggest the following uni- 
fied picture: The GRB central engine may have a range 
of magnetization parameter uq. 1. For low oo bursts, 
dissipation can occur at small radii, so that a bright 
photospheric emission component (such as the case of 
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GRB 090902B) is detected. Such bursts are usually ac- 
companied by a high energy component due to upscat- 
tering of the thermal photons (and probably also syn- 
chrotron self-Compton, Pe'er et al. 2012). 2. For inter- 
mediate (To bursts, the photosphere component is weaker 
but still detectable. Examples of this category include 
GRB 110721A (Axelsson et al. 2012) and GRB 100724B 
(Guiriec et al. 2011). The Band component of these 
bursts are formed in the large radii via internal shocks 
(IS) or internal collision-induced magnetic reconnection 
and turbulence (ICMART). 3. Finally, if o~q is large 
enough, both the photosphere and IS components are 
suppressed. The Band component forms at even larger 
radii via the ICMART process (Zhang & Yan 2011). 

Finally, very rapid hard-to-soft E p evolution observed 
in GRB 110721A (Axelsson et al. 2012 and this work) 
and many other GRBs (Lu et al. 2010, 2012) challenges 
existing models. Such a rapid evolution is not expected 
in the internal shock model. For the photosphere model, 
some moderate hard-to-soft evolution may be expected 
due to the intial growth of optical depth in a fireball (W. 
Deng & B. Zhang 2012, in preparation), but never an ex- 
treme evolution like GRB 110721A. A possible interpre- 
tion may be made within the framework of the ICMART 
model. According to this model (Zhang & Yan 2011), a 
in the emission region rapidly decreases during each IC- 
MART event, since the magnetic energy is continuously 
dissipated. If a good fraction of local magnetic dissipa- 
tion energy is deposited to electrons, the typical electron 
Lorentz factor j e would have a er-dependence, so that E p 
decreases with time as a reduces. Also relativistic turbu- 
lent reconnection would lead to locally Doppler-boosted 
mini-jets, whose Lorentz factors would also depend on 
the local a value (Zhang & Zhang 2012). The time de- 
pendent Doppler boosts for mini-jets would enhance the 
hard-to-soft E p evolution. 
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